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The hydrodynamic equation for the spatial and temporal evolutions of the electron temperature 
T e in the breakdown of the quantum Hall effect at even-integer filling factors in a uniform current 
density j x is derived from the Boltzmann-type equation, which takes into account electron- 
electron and electron-phonon scattering processes. The derived equation has a drift term, which 
is proportional to j x and the first spatial derivative of T e . When applied to the spatial evolution 
of T e in a sample with an abrupt change of the width along the current direction, the equation 
gives a distinct dependence on the current direction as well as a critical relaxation, in agreement 
with results of recent experiments. 
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The diagonal conductivity a xx is extremely small in 
the quantum Hall effectocP (QHE) at small currents. 
When the current is increased up to a critical value, a xx 
increases by several orders of magnitude within a—nar- 
row range of current, and the QHE breaks down-'C'?' 
(reviews are given in refs. |^, @, ||). The hysteresisdu'EP 
with the. f] h "i n gf pf the currenLaa well as the variety 
of spatiaoE3'li3Eif and temporaJu'EP behaviors observed 
in the breakdown of the QHE can be understood to be 
caused by a nonequilibrium phase transition. One of the 
most promising theories for describing the nnnecuiilib- 



rium phase transition is a hot-electron theoryfiolIEP in 
which the change of a xx is ascribed to the change of the 
electron temperature T e . In this research, the spatial 
evolution of <r xx in the current direction is studied on 
the basis of the hot-electron theory. 

The spatk^vplutions of a xx have been studied exper- 
imentallyoli3-liJ't3> in a system with an abrupt change of 
the sample width at x = (the x axis is taken along the 
current), which is designed to produce the change of the 
current density j = (j x ,0) at x — 0. When the sample 
is designed so that \j x (x < 0)| < j c i < j C 2 < \jx(x > 0)|, 
with j c i and jci the critical current densities (the system 
is bistable at j c \ < \j x \ < j C 2)i the relaxation in x > 
of a xx to the larger stable value and, the development of 
the breakdown have been observed.BliJ'LLS' On the other 
hand, when \j x (x < 0)| > j c2 > j c i > \jx(x > 0)|, the 
relaxation in x > of a xx to the smaller stable .value 
and the recovery of the QHE have been observedO A 
unique feature in this nonequilibrium phase transition is 
that the length scale for the relaxation of a xx depends 
strongly on the current direction: it is much larger for 
the electron flow in the +x direction. It is only for this 
direction of the electron flow that the relaxation length 
in x >, XL diverges as j x (x > 0) approaches a critical 
value. Ej'liilliiJ According to the hot-electron theory, this 
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current-direction dependence can be understood_by. the 
drift motion of hot electrons along the current.oO In 
order to understand the experiments and other spatial 
and temporal evolutions in a unified manner, the gen- 
eral theoretical framework is highly desirable. 

In this study, we derive the hydrodynamic equation 
describing the spatial and temporal evolution of T e from 
the Boltzmann-type equation on the basis of thejiiicro- 
scopic scattering processes proposed previously.E^P We 
restrict our derivation to even-integer filling factors 2N 
(N = 1, 2, • • •), for simplicity, and assume a uniform cur- 
rent distribution. We apply the obtained hydrodynamic 
equation to calculate the spatial relaxation of T e in the 
current direction. Our calculation reproduces both the 
current direction dependence and the critical behavior of 
the relaxation length. 

From a microscopic point of view, the potential in 
the plane of the two-dimensional electron system (2DES) 
fluctuates due to the presence of the layer with ionized 
donorsjLZP and the average distance between a potential 
hill and a neighboring valley, i v h, is Z v h ~ 0.1 /im 3> I 
with I the magnetic length. The screened potential 
has an amplitude, papproximately equal to the Landau 
level separation,13Eif and electrons (holes) populate the 
N + 1th (Nth) Landau level, even with a filling factor 
of 2N. Electron orbits in such potential consist of lo- 
calized closed orbits around a potential hill (valley) and 
extended orbits at the center of the broadened Landau 
level. Wave packets constructed from extended states 
perform drift motions perpendicular to the electric field. 
In addition, an electron hops between a large hill orbit 
and a large valley orbit near an extended orbit. Such 
a hopping motion is given by an electron-electron scat- 
tering, in which one electron hops between a hill and a 
valley and, at the same time, another is excited or re- 
laxed within a hill or a valley. Hopping motions along 
the electric field excite electrons, resulting in the energy 
gain of the 2DES. Electron-phonon scatterings, on the 
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other hand, relax the 2DES, resulting in the energy loss. 
In our previous study,t9 the energy gain and loss have 
been calculated for the higher T e branch and the van- 
ishing lattice temperature in homogeneous states. The 
theory is generalized in this study to include spatial vari- 
ations of T e in a macroscopic scale (much larger than l v ^) 
as well as an effect due to the drift motion of electrons 
in extended orbits (this effect is absent in homogeneous 
states). 

In constructing the Boltzmann-type equation, the 
2DES is divided into regions separated by an equipoten- 
tial line of the fluctuating potential so that each region 
contains one hill or one valley. The distribution of elec- 
trons is described by the occupation probability f a ia(t) 
of the ath state in the ith region at position in level a, 
where a is a combination of the Landau index and the 
spin and takes one of N f , N |, N+ 1 |, or N+ 1 j. The 
Boltzmann-type equation is then given to be 

dfaia ( $faia \ , / df ( 



Of 



V dt 



V dt 



(1) 



/drift \ /coll 

consisting of the drift term due to extended states and 
the collision term due to electron-phonon and electron- 
electron interactions. In electron-phonon scatterings, 
we have only included transitions aia —> aib within 
the same region (i) within the same Landau level (a), 
since their rates are much larger, than others due to 
larger overlaps of wave functions.EEP In electron-electron 
scatterings, we have only included transitions aia — > 
ajb, f3kc — ► pid within the same Landau level (a, 0). 
In the following we consider two classes of dominant 
electron-electron scatterings: (el) i = j and k = I, (e2) 
i = j, and k and I are the nearest neighbors. 

The local equilibrium is assumed in the scale of Z v h, 
since the inelastic scattering length due to electron- 
electron scatterings is estimated to be ~ I <C Z v h- Intro- 
ducing the electron temperature T e (ri, t) and the chem- 
ical potential p(ri,t) of the ith region, the occupation 
probability is given by 



where f(e\fi,T e ) = l/{exp[(e - /x)/fcsT e ] + 1}, and e a i a 
is the energy of state aia. 

A state in the local equilibrium is specified by the elec- 
tron temperature T e (r, t), the chemical potential p{r, t), 
and the electric field E(r,t) (these are assumed to be 
slowly varying in the scale of Z v i,)- The hydrodynamic 
equations for these variables are coupled in general. Here 
we assume that the current density is uniform, neglecting 
deviations near the sample edges. Although the current 
distribution in the Hall bar geometry is not uniform in 
the QHE,E3) it is suggested to be uniform in the tpeak- 
down region by the observations in many samplesEU'EHP 
that the critical current is proportional to the sample 
width. E is then given by E — (p xx j x , p yx j x ), with a 
constant current density j — (J x , 0) . Since p xx is smaller 
than p yx by, at least, one order of magnitude and p yx 
is quantized even in the breakdown regime, E can be 
approximated by E = (0, — j x h/2e 2 N). The chemical 
potential is then at the middle point between the iVth 



and the TV + 1th Landau levels and V/i = eE. Owing to 
the assumed uniformity of the current, the equation for 
T e (x,t) is decoupled from those for E and p. 

To obtain the equation for T e (x,t), we introduce the 
energy density U(x,t) defined by 

1 



U{x,t) 



(3) 



with pi = p(ri). An average is taken over regions within 
area S l^ h ) in which the variation of T e is negligible. 
The time derivative of U(x,t) is given in terms of the 
electronic specific heat C e = dU/dT e as 



dU_ dTe 

dt - ° e dt 



(4) 



and is also given using eq. (Q) as 



dux 



(5) 



dU _ /dU\ fdUY /dU 

dt - W drift + W col i + Ween' V dt J coU 

The three collision terms are the electron-phonon term, 
and the electron-electron terms (el) and (e2), respec- 
tively. 

The drift term is given by 

'dU\ dj* x 

dt /drift 9X 1 

with the energy flux jfj = (j^ x , 0): 

3Ux = 27^c/exi n x> 



(0) 



(7) 



where j 



di 



eEy/h is the electron flow averaged 



within S (3> ^h) a ^ ^ ne filling factor one, / cx = 
l/[exp(huj c /2kBT e ) + 1] is the occupation probability of 
extended states, and uj c is the cyclotron frequency. The 
electron-phonon term is given by 



<9f7\ p 1 



coll 



■ib 



= -Pl, 



(8) 



(2) where W v mi 



■i is the transition rate and Pr is the energy 

ct%a — >ib 1 w 

loss per unit area per unit time. When \e a i a — £ a ib\ <C 
ksTe (this inequality holds approximately in the higher 
T e branch (T e > Xr\) since \e a i a — e a i \ ~ 0.1huj c and 
k B T cl = 0.324^,0), we have f aia (l - f mb ) « f aia (l - 
faia) = -k B T e f'(e ai a) and obtain 

Pl = C L f e , (9) 

with T e — kBT e lfwj c and 

-f'{£ m a)]{e ai a - e aib )Wl ia ^ lb . (10) 



The electron-electron term (el) is given, when dT e /dx <C 
T e /l vh , by 



(11) 



'W\ el _ d_ dTe 

dt /coll ~ dxK<> dx ' 

where kq is the thermal conductivity due to (el) pro- 
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cesses. When \e a ia 



K = 



W: 



8S 

a, /3.z,fc, a, b,c,c£ 



£aib\ < ksT e , kq is given by 

f(e ma )}{-f'(e pkc )}Ae 2 Ar 2 



aia — >ib 



(12) 

with As = s aia - e ai b, Ar = \r t -r k \, and W^^kd the 
transition rate. The electron-electron term (e2) is given 

by 



9t /coll 



dx 



(13) 



where Pg — {—e)E-j n is the energy gain per unit area 
per unit time, and j n and jy = 0^,0) are the electron 
flux and the energy flux due to hopping processes, respec- 
tively. When \e a ia - e a ib\ < ksT e and eEl v h < ksT e 
(the second inequality holds at E ~ B c i since E c \ = 
20V/cm (B = 5 T)E9> and eE cl l vh = 0.03?iw c ), we ob- 
tain 



Pg = C G f cx (l - / cx )£ 2 , 
where = eE£/%u> c (I the magnetic length) and 



(14) 



4£ 2 S 

Similarly we obtain 



E ^"["/'(^a)]. (15) 



9 h 

J [7a: 



-CG/ex(l - /ex) 



j? 2 ar P 



4T e dx 

Finally we obtain the equation for T e : 



(16) 



(9;F 



Pg-Pl, (17) 



where the total thermal conductivity k is 

k(t 6 ) = Ko + c G / ox (i-/ cx )^. (is) 

When a constant density of states is assumed, Cl, 
Cg-i and Ko are approximately independent of T e . Esti- 
mates of Cl and Cq are given in our previous paper .E£P 
In estimating Cl, the deformation potential scattering 
by acoustic phonons and the vanishing lattice temper- 
ature were assumed; in estimating Cq, the summation 
over k and / was limited to the case where either k or I 
coincides with i. Calculated values of the critical elec- 
tric field E cl oc {C L /C G ) 1/2 and the energy dissipation 
CxxE^ oc Cl were in agreement with the observed values 
in their orders of magnitude. Following the calculation 
of Cq, here, we calculate kq, in which we limit the sum- 
mation over k to the nearest neighbors of i since the 
contribution to kq of larger Ar decreases as Ar~ 3 . The 
estimation shows that kq is only 5 percent of the total 
thermal conductivity, n. The T e dependence of k is ne- 
glected in the following, since dnfdT e — at E = E c \ 
andT e = T c i. 

Gurevich and MintsE-J introduced an equation for T e 
in their study of spatial and temporal variations of T e 
in the bistable region. In their equation they included a 
term proportional to j x and dT e /dx, which was ascribed 
to a contribution from the Seebeck effect. Although the 
drift term derived above is similar in form to their term 



[dfex/dx = (df cx /dT e )(dT e /dx)], it has a different ori- 
gin. The Seebeck effect vanishes in the present case of 
even-integer filling factors because the density of states 
is symmetric with respect to the chemical potential. 
In stationary states, the equation for T e becomes 

M(fe) Sr = + fe - 7 ^ 2/cx(i - /ox) ' (i9) 

in a dimcnsionless unit where 



M(T e ) = nhuj c /k B C L l 2 , 

V(Te) = V0jxfcx(l - /ox)! 1 ; 

?/o = eE cl uj c /2irlC L , 



(20) 
(21) 



(22) 

x = x/l, j x = {-j x )/j c i = E y /E cl , and 7 = 2.23342. 
The estimation gives ijq = 50 at B = 5 T. By chang- 
ing the variable as t — —x, the above equation becomes 
the equation for "position" T e as a function of "time" t 
of a particle with mass M{T e ) in a position-dependent 
friction rj(T e ) and a potential V(T e ) defined by 

- dV/df e = f e - 7ix/ex(l ~ /ex), (23) 

and plotted in Fig. 1(a). We have chosen this direction of 
"time" because T e with increasing x approaches a stable 
value in a simple exponential form and is convenient for 
the initial point. 

The current direction dependence is immediately ob- 
tained from the fact that r\ cx j x . The motion of the 
fictitious particle is much slower for the positive friction 
(jx > 0) than for the negative friction [j x < 0), repro- 
ducing the observed dependence of T e (x) on the current 
direction. 

A critical behavior is expected to occur at \j x \ = 
id, jc2- We consider in this study the critical behav- 
ior at \j x \ — id, since j c % — oo in our simplified calcu- 
lation of Pg and Pl- At \j x \ — j c \, the homogeneous 
steady solutions are f e = and f e = f cl = 0.324. We 
consider below a relaxation from T e > T c \ [x = 0) to 
T e = (x = oo) when \j x \ < j cl (\j x \ < 1). 

First we show the critical behavior qualitatively us- 
ing the above equivalent mechanical system. The initial 
"position" is chosen to be Tj n i t = T e (ii n it) ^ fcl- At 
T e ~ finit, /ox ~ and f e = f in i t exp[-(£-a % i„it)/M 1 / 2 ]- 
Therefore the initial condition at x — Xi n n is T e = Tj n i t 
and df e /dt = f^t/M 1 ' 2 > 0, and the "particle" goes 
down the slope of the potential as shown in Fig. 1(a). 
When friction r/ is positive (j x > 0), the motion of the 
particle slows down at the gradual slope at T e ~ T c \ and 
it takes a longer time to pass through the slope as j x — > 1 
since dV/dT e — » (j x — > 1). This means the divergence 
of the relaxation length of T e (x) as j x — » 1. On the other 
hand, no such critical behavior occurs when j x < since 
the friction is negative. 

Next we calculate the length scale of the slow variation 
of T e {x). In this calculation, the second-derivative term 
can be neglected because the decay length Id of d 2 T e /dx 2 
is given by Id = IMj r\ when dV/dT e = and is estimated 
to be Id = 5 /jm at B — 5 T and T e — T c ±. The equation 
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is in good agreement with the numerical ones. The re- 
laxation length L s is then defined by 

L s /e = Tr Vo ] x {h h 2 )- 1 / 2 . (27) 

L s diverges as (1 — jx)~ x ^ 2 as shown in Fig. 1(c), and 
becomes ~ 10 ^m at j x = 0.97. _ . 

The observed diagonal resistivity /C^aO remains to 
take large values within distance L° xp of approximately 
50 pm (j x ~ 0.97 and B — 4.8 T) from the position of 
the hot-electron injection, and then undergoes a steep 
decrease. The calculated relaxation length, L s , of T e 
(corresponding to p xx ) is in agreement with the observed 
relaxation length, L° xp , in the order of magnitude and 
the critical behavior. The present theory under the as- 
sumption of T e r~j T c i is not applicable to the observed 
steep decrease in the lower- p xx (lower- T e ) region. 

In conclusion, we have derived the hydrodynamic 
equation for the electron temperature T e from the 
Boltzmann-type equation. We have applied the obtained 
equation to the spatial evolution of T e and have repro- 
duced two observed results: the current direction depen- 
dence and the critical behavior of the relaxation length 
for T e . 

The author would like to thank Y. Levinson, N. 
Tokuda, and Y. Asano for valuable discussions. 



Fig. 1. (a) "Potential" 
kBTe/hoJc, jx = 
of x = x/£ with 
solutions of cq. ( 



V(T e ) defined in eq. (|| (T e = 
(—jx)/jcl = Ey/Efi). (b) T e as a function 
n defined in eq. fl22h (thicker lines: numerical 
), thinner lines: approximate analytical solu- 

as 



tions, eq. (|26[)). (c) Relaxation length L 3 defined in eq. 
well as ho, h,2, and Tq defined in eq. (E5h. 



is integrated to give 
x — Xq 



= ~Jx 



dT — 

To h(T,j x ) 



(24) 



with h(f,j x ) = f 3 {e l ' 2f + er 1 ' 2 ? + 2l- 1 j 2 x f 2 . f = 
T e (io) is the temperature at which h(T, j x ) takes its min- 
imum. The numerical integration gives T e (x) plotted by 
thicker lines in Fig. 1(b), which shows that the length of 
the slow- variation region L s increases as j x — ► 1. To find 
an analytical formula for L s , we replace h(T,j x ) by its 
Taylor series with respect to T around its minimum up 
to the second order: 



h(T,jx) ~ h (j x ) 



with h (j x ), h 2 (jx), and T (j 
analytical integration gives 

(h \ 1/2 
f e -f = -{-^-\ tan 



h 2 x )[f-f (j x )] 2 , (25) 
plotted in Fig. 1(c). The 



{hph 2 ) 



1/2 



-(x - Xq) 



(26) 



which is plotted in Fig. 1(b) by thinner lines; the curve 
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